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(3J[)! Abstract. We discuss explicit solutions of the logistic model with variable parameters. Classical 

data on the sunflower seeds growth are revisited as a simple application of the logistic model with 
periodic coefficients. Some applications to related biological systems are briefly reviewed. 
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1. Introduction 



We consider explicit solutions of logistic equations with time-dependent parameters. Particular 
attention is paid to the case of periodic coefficients which is especially important for many biological 
problems due to a natural periodicity of the Earth rotations. Models of this type can be used to 
mimic a population's response to seasonal fluctuations in its environment or to mimic a population 
with several discrete life-cycle stages [5], [7J, [8] [22]. As an example, sunflower plant growth 
data from the classical paper (21] are reviewed and analyzed from the viewpoint of the day-night 
periodicity. We dedicate this paper to Professor Doron Zeilberger on his 60th birthday — his 
■ nontraditional way of thinking has motivated one of the authors (S. K. S.) to review this classical 
qq topic from a modern perspective. 
O 

o 

2. Solution of the Logistic Model with Time-Dependent Coefficients 
^ The nonlinear equation for the logistic, or Verhulst's, model [6], [7J, [8], [21], [25], [26] 

H ; 

dx , , / _ x 



it a(t)x { l -W)) (21) 

with time-dependent coefficients a (t) and (3 (t) can be reduced to the linear equation 

by the standard substitution z — 1/x. Integration of (12. 2p with the help of an integrating factor, 

jt (" exp (/ a {s) ds )) = if exp (/ a {s) ds ) ' (2 - 3) 
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results in the explicit solution of the logistic model (in terms of a weighted harmonic mean): 

(2.4) 

corresponding to the initial condition x (to) = xq. A traditional expression is given by 



1 =c - sl «w da ( t q( g ) c j,; aw ^ dg , 1 

x (t) \J t0 (s) x (t c 



xq exp I / a (s) ds 
x (*) = 3 7 °/ ?i x (2-5) 



ft a (s) ( f s \ 
1 + xq I exp I / a (v) dv ] ds 

J to P ( s > \J t J 



(see also [5], [6], [7], [8], [12], [13], [H], [23] and references therein). Here, we would like to consider 
the most general functions a (t) and /3 (t) such that all integrals exist and 

'"' ft f (s) ds = f (t) (2.6) 



dt Jt0 

by a corresponding version of the fundamental theorem of calculus (see, for example, [H], |18j). 
The well-known solution with constant parameters a and /3 can be written as 

' I + ~ I) e- a ^\ (2.7) 



x(t) \x(t ) 13 
or 

x(t ) = (28) 

where [3 = lim^oo x (t) is the carrying capacity. 

3. Logistic Models with Periodic Coefficients 

If variable coefficients a (t) and (t) in the logistic equation (12. ip are periodic functions of time, 
namely, 

a(t + T)=a(t), 0(t + T)=p(t), (3.1) 

an important question arises about the basic dynamics and long-term behavior of a biological system 
under consideration. One may approach this matter in the following mathematical setting. If the 
period is T = t\ — t , the solution (12 .4p reads 

lb 

^i)" a + ^)' (3 ' 2) 
where parameters are given in terms of integrals over the period: 

exp I / a (s) ds ) dt, b = exp I — / a (t) dt J . (3.3) 

>t PW \Jt J V J t J 

Repeating this process indefinitely, we arrive at the following discrete map 

lb i v 

a + — , (3.4) 



where by the definition x (t n ) = x n . This recurrence relation can be rewritten as follows 

x n+1 - x n = (1 - b)x n+1 M - Y~Zl Xn ) ' ( 3 ' 5 ) 
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which can be thought of as a difference analog of Verhulst's model (12 .ip with constant parameters. 
We refer to ( 13. 4 p as a discrete logistic map (with a slight change of parameters it is also called the 
Beverton-Holt difference equation jl], [6], |23j). 

But the first order nonhomogeneous difference equation of the form 

z n+1 = a + bz n (3.6) 

with constant coefficients a and b has the following explicit solutions [6], [TJ: 

z n = a^-^ + b n z , if 6^1, (3.7) 

and 

z n = an + z , if 6 = 1, (3.8) 

which can be verified by a direct substitution. Therefore the discrete logistic map solution has the 
form 

1 a ( 1 a 



x n 1 — b ^ \xq 1 — b J ^ ^ 

for any initial data x . If b < 1, this sequence converges to the limit, say := lim^oo x n , given 
by 

1 " (3.10) 



and our solution (13. 9p takes a more convenient form 



Then 

X n Xq 

which implies that x n > x,^ for all n, if xq > x^ and vice versa. Finally, from (13. lip one gets 

b n lnb, (3.13) 



d f 1 \ 1 dXji Xq Xqg j n 



— = —+(-- — )b n . (3.11) 



(3.12) 



dft V x yi J x d/xXi xqX qq 

which reveals monotonicity properties of this bounded sequence in general. Namely, our sequence 
is strictly increasing, when x < Xqo, and strictly decreasing otherwise. 

Traditionally, a special question which arises in the theory of equations with periodic coefficients 
is the existence of periodic solutions [2], [17] . In the case of the logistic model under consideration, 
when an explicit solution is available, one can only obtain a continuous periodic solution by letting 
x\ = xo in our equation (13. 2p . Thus 

— = a + —, (3.14) 

x x 

which corresponds to the special initial condition given by the expression ( I3.10p . According to 
(13. lip the simple stability property holds, if b < 1, solutions with all other initial data converge to 
this periodic solution, or an attractor, as time goes to infinity with precise monotonicity properties 
described above. By (13.111) 

1 1 n l >r, (3.i5) 



Xn Vn V^O 2/0 
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which reveals the contracting map property for two different initial conditions x ^ y (see also [7] 
and 0). 

4. Two-Stage Compound Logistic Models 

In general, splitting a time interval [t , t 2 ] as the union of two successive subintervals, say [t , t 2 ] = 
[to,ti] U [ti, ^2] , one can write 

^ — + ^° * — + ^ (4 1) 

x(ti) a; (t )' x(t 2 ) x(t x y 

where the corresponding coefficients are given by our general expressions (13. 3p . Then 

lb 

where the 'compound' coefficients are given by the following analog of the addition formula 

a = ai + a &i, b = b b\, (4.3) 
which can also be derived directly from ( 13. 3D . It is worth noting the following group property 

1? ^y) (0 Z ) = O'^y) ' r^xi) (0 Z ) = O'^y) (44) 

and 

1 a \ / 1 ai \ _ / 1 ai + a bi 

6 J V 6 i / V b o b i 
(Another applications of the group of upper triangular matrices can be found in Refs. [IE] , [19] and 
[21].) 

With a slight change of notations, consider two step functions on the interval [0,ti] : 

n u\ / a o, if < t < t a (+ \ S Po, if < t < t . , 

a(t) = \ ai , ift <t<t 1 ^^j/^ ift <t<t 1 (46) 

and then define them by 'periodicity' for all t > with the period T = t\. We shall refer to this 
case as a two-stage compound logistic model (see also [12] for a special case). Evaluation of the 
integrals in ( I3.3P gives 

1 ~ Tl ■ ^ ~ T °) Tl h (A 7\ 

6 = Ton, (4-7) 



(4.5) 



where by definition 
As a result 



Pi Po 

r = e - ao ' , Tl =e- aiitl - to) . (4.8) 

^ = iim = s J M \7+« \ (4 - 9) 

Po + Pin - {Po + Pi T o) n 



n— >oo 



provided t ti < 1. 

Two special cases are of particular interest. If «i = 0, then t\ — 1 and 

1_r " 6 = r = e- Qoto , x oo = /3 . (4.10) 



A) 

This will be useful in section 8. Another limit, ( 5 1 — >■ oo, corresponds to a resetting mechanism 
discussed in [15] (see also [23]). 



the logistic model 
5. w-Stage Compound Logistic Model 



5 



Results from the previous section can be easily generalized in the following fashion. Condider a 
partition t = s < $i < ■■■ < s n = t\ of a time interval [to, t±] = U^Tq 1 [sj, Si+jJ . Then with the help 
ofgSD-gS}: 

1 b 



where 



or in the matrix form 

n—l 



x(s n ) ~ 1 x(s y 

n—l n—l n—l 

J> II 6 *> b = U b ^ ( 5 - 2 ) 

i=0 k=i+l i=0 



i=0 



nn-l 
j=n 



by induction. It is worth noting that this multiplication property holds for any numbers {a{\ 



and {bi}™ = Q , thus generating solution of the following difference equation 



n-1 
i=0 



' a l + ^, (5.4) 



which is obtained here with the help of the group property of upper triangular matrices (see also 
[7j and [23] )• In this paper, we take = x (si) and use definitions f)3.3p . namely, 

f s,+1 a (t) ( /"* \ / f Si+1 \ 

a>i = bi / — — — exp I / a (s) ds J c?i, = exp I — / a (t) dt ) (5.5) 



/3(t) 

for the corresponding subintervals [s^ Si+i] with i = 0, 1, n — l. 

For step functions a (t) = a, and /3 (i) = /3 i on subintervals Si < t < s i+ i one gets 

1 - r. 



and 



^ = Ti = e - a <( fl *+i-*) (i = 0,1,..., n-1) (5.6) 



n—l n—l n—l 

i=0 P % k=i+l i=0 



Further details are left to the reader. 



6. An Extension 



The following generalization of the logistic model 

dx ( x e 



dx f 

a (t) x ( 1 — ) , = constant (6.1) 
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has been considered in Refs. [T], [3] and [6]. This equation is integrable by the substitution z — l/x 6 
and the corresponding extension of (12. 4p takes the form 



> ./£«(,) ds ( n f a{s) 0f: o a(v) dv . ' 



with x — > x e and a — > 9a. The case of periodic coefficients is analyzed in a similar fashion with the 
help of this substitution. Further details are left to the reader. 



7. Logistic Models and Riccati's Equation 



Considering a more general model 

l = o w*( 1 -^o) +7(t) " (71) 

where a (t) , (3 (t) and 7 (t) are suitable real valued functions of time only, one obtains the so-called 
Riccati equation, which can also be reduced to a linear second order differential equation by a 
certain substitution (see, for example, Refs. [9] and [10] for more details). The case of periodic 
coefficients will be discussed elsewhere. 



8. An Application: Logistic Growth of Sunflower Seeds Revisited 

Logistic models are well-known in auto catalysis, population dynamics, mathematical epidemiol- 
ogy and theory of fishing (see, for example, [1], [6], [7], [8], [20], [23] and references therein). Among 
other biological applications, we utilize the periodic logistic models, described above in detail, to 
sunflower plant growth data from the classical paper [21], where the sunflower (helianthus) was 
chosen because of the fact that it grows without producing branches and it was thought that mea- 
surements of height and weight represented the growth of the entire organism with a fair degree of 
accuracy. Original experimental results are presented in the following table: 



Sunflower height versus growing days by Reed and Holland (1919) 



days n 


observed mean height x n (cm) 


constant i^(this paper) 


In ( Xn ) 





10.00 


0.040593 


-3.1966 


7 


17.93 ±0.14 


0.04119 


-2.5798 


14 


36.36 ± 0.43 


0.03852 


-1.7917 


21 


67.76 ± 0.78 


0.033353 


-1.0137 


28 


98.10 ± 1.38 


0.032672 


-0.46643 


35 


131.00 ± 1.73 


0.032005 


0.058956 


42 


169.50 ±2.21 


0.038430 


0.6902 


49 


205.50 ±2.92 


0.042069 


1.4336 


56 


228.30 ±3.41 


0.043129 


2.1649 


63 


247.10 ±3.80 


0.052904 


3.5083 


70 


250.50 ± 3.76 


0.050188 


4.1372 


77 


253.80 ±3.99 


0.059799 


5.8932 


84 


254.50 ± 3.89 


not available 


not available 
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The following solution of the continuous Verhulst model, see (I2.7p - fl2.8p . 



with = 0.0421 and t\ = 34.2 was originally suggested in order to interpolate these data [21] . Here, 
the earring capacity was chosen to be /3 = 254.5 and the time at which the growth has run half way 
to equilibrium; that is the time at which x = (3/2 and the maxinum growth occurs; was estimated 
as t\ = 34.2 days. It is worth noting that the average value of the constant K, reevaluated here 
from our corrected version of the third column (without n = 0), is given by K = 0.042205, which 
is slightly different from the value K = 0.0421 obtained in the original paper [21] (this average 
becomes K = 0.042071 « 0.0421 if our n = data is included). 

A simple assumption of the day-night periodicity of the plant growth results in our solution (13. lip 
of the discrete Verhulst map, which can be rewritten for the numerical analysis purposes in a more 
convenient form: 

In ( Xn ) =\n( Xl ) +(n-7)\n({) . (8.2) 



Then we match the data available from the table on the seven day scale only from x 7 = 17.93 till 
£§4 = 254.50 ~ Xqc,. Corresponding data from the last column have been first analyzed with the 
help of a least square linear curve fitting on the entire available time interval [7, 77]. As a result 

In ( — ) = -3.612328764 + 0.112057018181818180™ (8.3) 

V254.50 -x n y v ; 

= 0.112057018181818180 (n - 32. 237) , 

which gives the following modified values of parameters t\ = 32. 237 and K = 0.112057018181818180 
/ In 10 = 0.048666 in the original formula (18. ip . Unfortunately, this formula does not match observed 
data better then the original expression in Ref. [2T] (even when one uses a quadratic least square 
approximation for the last column; two (or several)-stage logistic models seem more approapriate). 
An elementary data analysis of the last column shows that the best linear approximation occurs 
during the time interval [7, 49] , when by the least square method: 

In f — ] = -3.106504286 + 0.0922278367346938688n (8.4) 

V254.50 -x n / 

= 0.0922278367346938688 (n - 33.68293561) 

with the best values of parameters t x = 33.68293561 and K = 0.04005404057 in the formula (]8TT) . 
Original data are approximated well by sigmoid shape function; the revised results are presented in 
the following graph (see Figire 1). 

Finally, we would like to apply our periodic two-stage (day-night) compound logistic model to 
helianthus growth. One may speculate that the sunflower plant was growing about fourteen hours 
during an average California day from May to August of 1919 and was not growing at all during 
nights. Then by fl4.10p : 

7 



exp ( -j^®o 



bT\ b = (8.5) 



If the (sun)light would be provided for twenty four hours every day, say in a greenhouse, expression 
(18.1 p should be modified as follows 

l0 « 254^ = ^ <"-*>> < 8 ' 6 > 
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with K = (12/7) K = 0.068664 and t = (7/12) t x = 19.648. The corresponding graph is also 
presented on Figure 1. 

Height Comparison with Best Estimate Fit Parameters 
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Figure 1 . Comparison of observed and calculated values for the mean hight of helianthus. 
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